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1.  Introduction 


This  paper  will  detail  the  research  done  on  the  non-linear  refractive  index  of 
various  optical  materials,  specifically,  Bismuth  Silicon  Oxide  (BSO),  and  a  Ce-doped  Ba2- 
xSrxKi-yNayNbsOis  (KNSBN;Ce)  crystal.  The  areas  presented  are  the  use  of  the  Z- 
scan  technique  to  measure  the  non-linear  refractive  index  of  the  crystals,  use  of  the 
samples  as  optical  limiters,  investigation  into  the  self  phase  modulation  inherent  to  the 
samples,  and  experimental  results  for  each  of  the  areas.  Recently,  the  self  focusing  in  new 
optical  materials  such  as  KNSBN:Ce  and  bacteriorhodopsin  (BR)  under  low  CW  HeNe 
illumination  were  reported! In  this  paper,  the  properties  of  the  BSO  and  KNSBN:Ce 
crystals  were  investigated  using  a  high  power  Argon-Ion  laser. 

The  major  property  of  interest  is  that  of  the  non-linear  refractive  index  of  the 
crystals.  Because  of  this  nonlinear  refractive  index,  a  Gaussian  beam  propagating  through 
the  crystals  induces  a  lens-like  refractive  index  profile.  This  refractive  index  change  in 
turn  modifies  the  beam  propagation.  It  is  important  to  understand  this  lensing  effect  in 
order  to  better  use  the  material  in  optical  system  applications. 

The  Z-scan  technique  offers  a  non-complicated  yet  precise  method  of  measuring 
the  effective  non-linear  refractive  index  coefficient  of  the  samples.  All  that  is  needed  is 
the  laser,  a  focusing  lens,  an  aperture,  and  a  power  meter.  This  method  is  based  on  the 
self-focusing  or  defocusing  phenomena  that  occurs  as  the  sample  is  translated  through  the 
focus  of  the  lens  in  the  direction  of  propagation,  or  the  Z  direction. 

After  the  measured  and  derived  results  are  presented,  the  use  of  the  samples  as  an 
optical  limiter  is  investigated.  The  optical  configuration  is  the  same  except  that  the 
crystal  is  stationary  while  the  power  is  increased.  An  aperture  is  used  in  conjunction 
with  the  self-focusing  or  defocusing  effect  of  the  crystal  to  limit  the  amount  of  power 
transmitted  down  the  optical  path. 

Finally  the  effect  of  self  phase  modulation  on  the  samples  is  investigated.  Self 
phase  modulation  occurs  along  with  self  focusing/defocusing  and  has  an  important 
influence  on  the  optical  limiting  behavior  of  the  crystals.  A  comparison  between 
experimental  and  measured  data  is  also  presented. 

2.  Z-Scan  Theory 

The  Z-scan  is  a  sensitive  and  convenient  method  for  measuring  the  light  induced 
changes  in  the  effective  nonlinear  refractive  index  and  absorption  coefficient  of  optical 
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materialsl^  The  Z-scan  is  based  on  the  transformation  of  the  phase  distortion 
associated  with  the  self-lensing  into  an  amplitude  distortions  during  the  beam 
propagation. 

The  experimental  arrangement  is  shown  in  Fig  1.  A  lens  is  used  to  focus  the 
illuminating  laser  as  the  crystal  is  translated  along  the  optical  path  (the  Z-scan).  A 
focusing  lens  is  used  to  suppress  the  beam  farming  effect  of  photorefractive  materials, 
which  is  proportional  to  the  beam's  lateral  dimension,  and  to  make  the  self  focusing  the 
dominant  effect.  In  particular,  ICNSBN  has  a  strong  fanning  effect  even  under  low 
illumination  intensitiesl^l  A  limiting  aperture  is  placed  past  the  focus  of  the  lens  and  a 
detector  is  used  to  measure  the  output  intensity.  The  aperture  is  set  to  let  only  half  the 
light  intensity  pass  when  no  crystal  is  in  the  path.  The  aperture's  intensity 
transmittance,  s ,  is  defined  as  the  ratio  of  output  power  to  input  laser  beam  power.  The 
crystal  is  placed  directly  behind  the  lens  and  translated  through  focus.  The  defocusing 
(focusing)  of  the  crystal  will  allow  more  (or  less)  of  the  light  to  pass  the  aperture  by 
creating  a  longer  (or  shorter)  effective  focal  length  of  the  lens  and  crystal  system.  The 
beam  diameter  at  the  aperture  is  compressed  or  expanded  with  this  changing  effective 
focal  length.  Measurements  of  the  intensity  at  the  detector  with  respect  to  crystal 
position  give  the  position  dispersion  curves  which  can  be  used  to  find  the  sign  and 
magnitude  of  the  focusing.  A  pre-focal  transmittance  minimum  (valley),  followed  by  a 
post-focal  transmittance  maximum  (peak)  indicates  a  positive  lensing.  Fig.  1  shows  the 
case  of  positive  dispersion,  with  less  then  more  light  transmitted  by  the  aperture  as  the 
sample  is  translated  through  focus. 

It  is  believed  that  the  focusing  effect  in  KNSBN  and  BSO  is  due  mainly  to  thermal 
effects  of  the  incident  radiation  on  the  index  of  refraction  of  the  crystaU'l  The  intensity 
dependence  of  the  refractive  index  can  be  written  as  n  =  +  An  =  n„  -I-  n2l^  where  n„  is 

the  index  of  the  crystal.  An  is  the  light  induced  index  change,  and  is  the  effective  non¬ 
linear  refractive  index.  The  magnitude  of  nj  can  be  determined  froml^l 
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where  is  the  difference  between  the  peak  and  valley  of  the  position  dispersion 

curve,  oc  is  the  linear  absorption  coefficient,  h  is  the  thickness  of  the  sample,  f  is  an 
empirical  constant  defined  as  /  =  0.406(1  -  where  5  is  the  aperture  linear 
transmittance,  and  is  the  intensity  at  the  Gaussian  beam  waist  inside  the  sample 
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200  mm  Aperture 

(c) 

Fig.  1  Experimental  setup  for  Z-scan  measurement. 


The  Z-scan  can  be  calculatedf’^  by  knowing  that  the  power  transmitted  through 
the  aperture  of  radius  is  given  by 


PAz)  =  P, 


2r 


l-exp(- 

o). 


t) 


(2) 


where  is  the  beam  radius  at  the  aperture  and  is  the  power  at  the  aperture.  If  the 
power  fluctuates  with  time,  the  transmitted  power  is  given  byt^i 


jPTiz,t)dt 

nz)=^ -  (3) 

sjPa(l)dl 

0 

where  s  is  the  aperture  transmittance.  If  temporal  variations  can  be  ignored,  the  above 
equation  can  be  rewritten  ast^l 


nz)= 


2r 

l-exp( - f-) 

(o: 


(4) 


where  is  given  by  1^1 


2A%(D-x)x^ 

+ 

f  2A<I>o(Z)-x)^ 

col  ' 

[  aD(l  +  x^f  J 

1  a(l  +  xy  j 

(5) 


where  AOq  is  the  on  axis  nonlinear  phase  shift,  D=d/z^  and  x=  z/Zq  and  a  is  a  correction 
factor  which  a  suitable  choice  for  is  a  =  6.4(1  —  .  Fig.  2  shows  the  calculated  curve 

for  the  Z-scan  transmittance  compared  to  the  actual  data  for  KNSBN  at  488nm. 
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2.5 


z 

Fig.  2  Z-scan  for  KNSBN  at  488nm,  a  :  calculated,  b  :  measured. 


The  Z-scan  theory  presented  here  assumes  the  crystal  can  be  treated  as  a  thin 
sample  when  its  thickness  is  less  than  the  diffraction  length  of  the  focusing  lens.  The 
diffraction  length,  z^i,  is: 


Taking  the  data  in  the  experiment  using  BSO,  (O  is  32|j,m,  n  is  2.6,  and  ^  is  5 14.5nm,  the 
value  of  z^i  is  approximately  14mm.  The  BSO  crystal  thickness  is  4mm  so  it  can  indeed 

be  treated  as  a  thin  sample. 

When  the  length  of  the  crystal  is  on  the  order  of,  or  greater  than  z^i,  the  Z-scan  in 
Fig.  3  can  be  expectedl^-^f  The  peaks  and  valleys  of  the  Z-scan  trace  have  become  more 
pronounced  and  further  separated.  A  conceptual  understanding  can  be  obtained  if  the 
flattened  transmittance  curve  about  Zg  is  considered  first.  When  the  thick  sample  is 
centered  at  Zg ,  the  crystal  will  create  large  local  phase  distortions  in  the  beam,  but  the 

prefocal  distortion  will  be  balanced  out  by  the  postfocal  distortion.  It  can  be  thought  of 
as  two  lenses  centered  about  Zg  with  the  second  lens  canceling  the  effect  of  the  first  lens. 

Therefore,  the  far  field  pattern  will  be  relatively  unchanged  giving  a  flat  transmission  line. 

As  for  the  peak  and  valley,  the  maximum  and  minimum  occur  when  the  focus  is  at 
the  surfaces  of  the  sample.  Take  the  case  of  positive  n2-  When  the  focus  is  at  the  back 
face  of  the  sample,  all  of  the  lensing  is  acting  to  decrease  the  system  focal  length  leading  to 
a  minimum  in  the  transmission  curve.  As  the  sample  begins  to  move  through  focus,  the 
competing  effects  of  the  prefocal  and  postfocal  phase  distortions  start  to  cancel  each 
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other,  leading  to  the  flattened  curve.  The  opposite  happens  after  focus  when  the  focal 
point  is  at  the  front  surface  of  the  crystal  and  all  the  lensing  tends  to  focus  the  beam, 
allowing  more  light  to  pass  the  aperture  leading  to  a  maximum  in  transmittance. 

For  the  thin  sample,  the  peak  and  valley  are  separated  by  approximately  1.7  z^i- 

For  the  thick  sample,  since  the  thickness  is  now  on  the  order  of,  or  greater  than,  the 
diffraction  length  the  thickness  dominates  and  begins  to  push  the  peak  and  valley 

apart.  The  peak  and  valley  occur  when  the  focus  is  at  either  the  front  or  back  surface. 

The  solid  line  in  Fig.  3  is  for  the  case  of  positive  n^,  and  the  dashed  line  is  for  negative  n^. 


Fig.  3  Position  dispersion  curves  for  a  Z-scan  of  a  thick  sample. 


3.  Z-scaii  experimeHtal  results 

Position  dispersion  curves  for  a  sample  of  BSO  are  shown  below  in  Fig.  4.  These 
position  dispersion  curves  can  be  used  to  obtain  a  value  of  Uj  BSO.  Each  of  the 
curves  correspond  to  a  different  input  power  level  which  is  the  reason  for  the  different 
amplitudes  of  the  curves.  Note  that  in  Eq.  1,  ^2  is  related  to  the  difference  in  peak  and 
valley  divided  by  the  input  intensity  so  that  each  position  dispersion  curve  should  yield 
similar  values  of  From  these  curves,  the  average  was  found  to  be  1.47  x  lO'^cm^AV . 

This  is  comparable  with  the  Uj  =  1.9  x  lO'^cm^/W  reported  in  Ref.  7. 
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z  (  mm  ) 

Fig.  4  Position  dispersion  for  the  BSO  crystal  for  powers:  curve  a,  0.1  mW;  b, 
1  mW;  c,  5  mW  and  d,  10  mW. 


For  the  KNSBN:Fe  crystal,  position  dispersion  curves  are  presented  in  the  Fig.  5. 
The  figures  are  for  5 14.5nm,  488nm,  and  476nm.  It  is  noted  that  the  magnitude  of  the 
curves  varies  with  the  polarization  of  the  incident  illumination.  The  polarizations  were 
parallel  and  perpendicular  to  the  C-axis  of  the  crystal.  This  can  be  explained  by  a 
difference  in  the  absorption  coefficients  of  KNSBN  for  orthogonal  polarizations!' I  The 
absorption  affects  the  change  in  the  index  of  refraction  leading  to  a  larger  or  slightly 
smaller  position  dispersion  curves.  In  other  words,  the  less  absorption  the  crystal  has, 
the  smaller  the  change  in  the  index  of  refraction  An,  leading  to  a  smaller  variation  in  the 
peak  and  valley  of  the  position  dispersion  curve.  For  the  opposite  case  of  higher 
absorption,  the  position  dispersion  curve  would  naturally  exhibit  a  greater  variation  in  the 
peak  to  valley.  Table  1  below  gives  the  values  of  the  absorption  for  the  two  polarizations 
at  the  wavelengths  used.  The  for  KNSBN  is  positive.  Note  that  this  is  different  from 
KNb03  where  the  light  induced  lensing  is  negativeP  '®!.  Also  a  definite  difference  can  be 
seen  in  the  absorption  coefficients  for  the  different  polarizations. 
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n2(cm2/W)  1.06x10-7  i. 31x10-7  0.77x10-7  1.08x10*7  0.79x10-7  1.03x10-7 


Table  1  Results  of  KNSBN  at  514.5nm,  488nm,  and  476nm. 


4.  Optical  Limiting 


An  optical  limiter  is  a  device  which  produces  an  output  power  that  is  decreasingly 
proportional  to  the  input  power  when  the  input  is  beyond  a  certain  threshold  value.  It 
has  many  applications  in  laser  damage  protection.  The  self-focusing  nature  of  KNSBN 
and  BSO  can  be  used  to  limit  the  optical  intensity  transmitted  through  a  systemt^.^.s].  By 
using  a  lens  and  an  aperture  as  shown  in  Fig.  1,  and  placing  the  crystal  in  the  valley  of  the 
position  dispersion  curve,  the  intensity  transmitted  is  not  allowed  to  increase  linearly  but 
is  limitedt^l.  As  the  laser  power  is  increased,  the  intensity  measured  before  the  aperture  is 
taken  to  be  the  input  power  and  the  intensity  past  the  aperture  is  the  output  power.  The 
plot  of  output  vs  input  power  shows  the  nature  of  the  crystal  optical  limiting.  The 
optical  limiting  experiment  was  repeated  for  various  apertures;  it  was  found  that  the  value 
of  5  around  0.5  produces  the  best  optical  limiting  curve,  i.e.  the  output  is  very  flat.  This 
experimental  data  supports  selection  of  5  to  be  0.5  in  the  Z-scan  experiment.  Fig.  6 
shows  various  optical  limiting  curves  for  the  BSO  crystal. 

For  KNSBN,  the  limiting  was  measured  at  three  Argon  wavelengths,  5 14nm, 
488nm,  and  476nm.  It  can  be  seen  in  Figs.  7  and  8  that  the  optical  limiting  is  strongly 
wavelength  dependent.  This  is  most  likely  explained  by  an  absorption  wavelength 
dependence.  A  characteristic  dip  in  the  transmitted  intensity  is  evident  at  30mw.  This 
oscillation  can  perhaps  be  best  explained  by  self  phase  modulation  of  the  incident  beam. 

5.  Optical  Limiting  experimental  results 


0  2  4  6  8  10 

Input  Power  in  mW 

Fig.  6  Behavior  of  BSO  optical  limiting  for  various  aperture's  linear  intensity  transmittance: 
curves  a,  b,  c  and  d  are  for  ,?  values  of  0.2,  0.5,  0.66  and  0.8  respectively. 
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0  5  10  15  20  25  30  35 

Input  (mW) 

Fig.  7  Optical  limiting  for  KNSBN  for  514nm,  488nm  and  476nm  up  to  first  dip. 


0.1  1  10  100  1000 


Fig.  8  Optical  limiting  for  KNSBN  for  514nm,  488nm  and  476nm  at  higher  powers  and  simulation  data. 

6.  Self  Phase  Modulation 

Self  phase  modulation  (SPM)  occurs  in  conjunction  with  self  focusing.  SPM  is  an 
intensity  dependent  effect  and  is  related  to  the  non-linear  coefficient  rij.  As  the  beam 
propagates  through  the  crystal,  the  intensity  distribution  of  the  beam  affects  the  index  of 
the  crystal  which  in  turn  affects  the  propagation  of  the  beam.  The  change  in  the  index  of 
refraction  not  only  has  a  focusing  effect  on  the  beam  but  also  causes  a  phase  change  in  the 
beam  as  well.  The  crystal  will  appear  to  have  a  different  thickness  corresponding  to  the 

1  1 


change  in  the  index  of  refraction.  In  the  context  of  self  focusing,  the  Gaussian  beam 
profile  of  the  incident  radiation  leads  to  a  varying  phase  change  across  the  beam  profile 
given  by 


=  (7) 

where  lo/r)  is  the  Gaussian  illumination  profile  and  is  the  effective  optical  thickness 
defined  as 


where  a  is  the  effective  absorption  coefficient. 

As  the  intensity  increases  to  the  point  that  the  phase  change  exceeds  2n,  the  far- 
field  pattern  begins  to  exhibit  a  concentric  ring  pattern.  More  rings  become  apparent  as 
the  phase  continues  through  multiple  27t  phase  shifts.  The  total  number  of  rings  N  can  be 
estimated  by  a  following  relationship, 


yy _  max 

2k 


n2iJe 


(9) 


The  far  field  projection  of  this  complex  intensity  distribution  can  be  calculated  by 
assuming  the  input  intensity  is  Gaussian  and  by  using  the  phase  shift  expressed  in  Eq.  7. 
The  field  amplitude  at  the  exit  surface  of  the  crystal  can  be  expressed  asl"! 

E{r)  =  {I^{r)f^e\^{i[k,z  +  ^{r)  +  (f)Q\}  (10) 

where  0(r)  is  the  phase  shift  from  Eq.  7  and  is  the  uniform  retardation  due  to  the 
of  the  crystal.  Simply  by  taking  the  Fourier  transform  of  the  above  equation,  the  far-field 
radiation  pattern  can  be  obtained. 

As  the  intensity  of  the  irradiation  increases,  the  phase  change  becomes  stronger 
which  varies  the  interference  pattern.  When  the  pattern  goes  through  a  minimum  at  the 
center  of  the  beam,  the  effective  transmittance  of  the  system  is  lessened.  The  oscillatory 
nature  of  the  phase  change  with  intensity  explains  the  oscillations  in  the  plots  of  the 
output  power. 
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A  simulation  of  the  above  theory  was  done  using  MatLab  and  the  results  were 
consistent  with  the  observed  experimental  results.  Appendix  B  contains  the  MatLab 
code  along  with  plots  of  the  intensity  profile  and  optical  limiting  curve.  Fig.  8  shown 
above,  contains  a  comparison  of  the  calculated  and  measured  optical  limiting  curves  for 
the  case  of  KNSBN  at  514nm. 

7.  Self  Phase  Modulation  Experimental  Results 

Fig.  9  shows  line  scan  plots  of  the  far-field  intensity  pattern  resulting  from  the 
self  phase  modulation  of  a  BSO  crystal.  Note  that  the  powers  used  in  this  figure  are 
much  higher  than  that  of  the  optical  limiting  curves  in  Fig.  6,  which  explains  why  no 
oscillations  in  the  optical  limiting  curves  are  apparent  Fig.  10  shows  the  predicted 
intensity  profiles  resulting  from  self  phase  modulation. 


Fig.  9  Transverse  intensity  profiles  of  transmitted  beam. 
Curves  a,  b,  c,  d  and  e  are  for  5,  53,  66,  136,  156  mW  respectively. 
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Fig.  10  Results  of  calculations  for  the  beam  profile  due  to  self  phase  modulation. 


The  results  in  Fig.  1 1  were  obtained  with  a  Sony  XC-75  black  and  white  camera 
and  a  Spiricon  beam  profiler  which  was  used  to  transfer  the  pictures  via  GPIB  to  an 
Apple  Mac  Ilex.  The  pictures  clearly  show  the  ring  structure  which  results  from  the 
phase  distortion  of  the  crystal.  These  results  were  from  a  BSO  crystal. 


Fig.  11  Far-field  ring  patterns  for  incident  intensities:  150mW,  302mW,  378mW,  and  514mW. 


8.  Conclusion 


The  experimental  method  for  the  measurement  of  the  non-linear  refractive  index  of 
BSO  and  KNSBN  has  been  presented  along  with  the  measured  results.  For  BSO  it  was 
found  that  nj  is  on  the  order  of  10"^  cm^AV.  This  result  is  consistent  with  previously 
reported  results  based  on  the  interferometric  method.  For  KNSBN  it  was  on  the  order  of 
10'^  cm^AV.  The  use  of  these  materials  as  optical  limiters  has  been  investigated  along 
with  the  best  configuration  for  their  use.  In  the  case  of  BSO,  it  was  shown  that  its 
performance  as  a  limiter  is  at  optimum  when  the  aperture's  linear  transmission  s  is  set  to 
0.5.  Both  crystals  pointed  out  that  for  higher  input  power  the  self  phase  modulation 
effect  could  not  be  ignored.  The  self-phase  effect  leads  to  a  concentric  ring  pattern  in  the 
far-field  when  the  maximum  phase  distortion  exceeds  2tz.  This  ring  pattern  has  a 
profound  effect  on  the  optical  limiting  curve.  Calculated  optical  limiting  curves  for  self 
phase  modulation  were  in  agreement  with  measured  data. 
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Appendix  A 


MacPhase  Macro  for  Calculating  Self  Phase  Modulation 


macro  SPM; 
var 

boolean  :  ok; 
longint:  i,  maxNum; 
real:  R,  rPhase.sum; 

Str255  :  gaussName,  lineName; 

begin 

DisposeUnClose; 

FFTSetting(TRUE, TRUE, TRUE); 

ok:=FrontData{gaussName,'data'); 

if(!ok) 

beep; 

exit; 

endif; 

maxNum:=50; 

lineName:='Line'; 

NewData(lineName,1  ,maxNum,0, false); 

MoveWindow(lineName,  400, 400, TRUE); 
for(i,1  ,maxNum) 

R:  =  EvalNumber(i,*,1 ); 
rPhase:=EvalNumber(i,*,0.25); 

SimpleMath(gaussName,*, 'None',  R,' Amplitude'); 
SimpleMath(gaussName,*, 'None', rPhase, 'Phase'); 
ok:=Exist('Farfield'); 
if(ok) 

DisposeData{'Farfield'); 

endif; 

FFT2D('Ampiltude','Phase',  128,1 28, 0,'Farfield','Farfieid  Phase'); 
SetROIOval{'Farfield', 59, 59,72,72); 
lntegrate{'Farfield',3,'asdf'); 

GetRindex(1  ,sum); 

PutDataNumber(lineName,1  ,i,sum); 

PlotData(lineName,'Line  Plot',TRUE); 

MoveWindow('Farfield',  400, 200, TRUE); 

PlotData('Farfieid', 'Color  Contour  Plot',TRUE); 

DisposeDataC  Amplitude'); 

DisposeData{'Phase'); 

DisposeData('Farfield  Phase'); 

endfor; 
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Appendix  B 


Matlab  for  SPM 


%  Plots  the  SPM  farfield  beam  pattern  and 
%  finds  the  power  in  the  central  half 
clear 

cig.hold  off 

load  argonX.load  argonY  %  loads  actual  data 

echo  on 

count=0; 


thehalf  =1253; 
x=(-64:. 05:63); 
alpha=4.0948; 

Lambda=514.5e-9; 
d=0.05225*1e-2; 
n2=1.06e-7; 
aa=(6.8e-4)^2; 
phaseF=(2*pi*d*n2)/(Lambda*aa)  %  as  in  Eq.  7 


O/ 

/o 


%  previously  determined  from  below 

%  absorption  KNSBN  in  cm^-1 
%  wavelength 
thickness  of  crystal 
%  from  Table  1 

%  spot  size  of  beam  in  crystal 


%  Guassian  profile 
sigma=1 ; 

f=sqrt(1/(2*pi*sigma^2))*exp(-(x-4).^2/(2*sigma^2)); 

norm=sum(sum(f)); 

f=f./norm;  %  make  integrated  area  equal  to  one 

plot(x,f) 

axis([-11  10  0  lei]) 
run=1 ; 

for  k  =  1:15:300,  %  power  in  mW 

k; 

kk=k*1e-3 
count=count+1 ; 
xx(count)=kk; 

phase=(phaseF*kk*f);phaser=max(max(phase))/pi 
U=((kk.*f).^0.5).*exp(i.*phase); 
UU=fftshift(fft(U));IUU=(UU.*conj(UU))/1; 
plot(x,IUU),hold  on 

if  run  ==  1,  %  Find  the  half  power  points 
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u2=iuu: 

[y,themax]=max(u2); 

y2=y/2; 

if  thehalf  ==  0, 
for  i=1  :length(u2), 
if  u2(i)>y2, 
u2(i)=0; 
end; 

if  i>length(u2)/2, 
u2(i)=0; 
end; 

[yy,thehalf]=max(u2);  %  Find  the  half  power  points 
end;  %  for  max 
end; 


run=0; 

diff=abs(themax-thehalf)/1 ; 

%  Find  the  power  contained  in  the  center 
b(count)=sum(IUU(themax-diff:themax+diff)); 

end; 

%pause 
hold  off 
%plot{b') 
bb=b./max(b); 

%  plot  calculated  data  with  actual  data 
plot((1e-3.*argonX),  argonY),  hold, piot(xx,bb), hold 


%  SPM  in  2-D 

%  July  28  1994  -  Self  Phase  Modulation 
clear 

s=32;step=.5; 
x=(-s:step;(s-step)); 
sigma=1 ; 

ax=0;ay=0;  %  to  offset  gaussian  from  center 

y=x'; 

X=ones(y)''x; 

Y=y*ones(x); 

R=sqrt((X-ax).^2+(Y-ay).^2)  ; 

f=sqrt(1/(2*pi*sigma^2))*exp(-(R).^2/(2*sigma^2)); 

for  i  =  1:1, 
num=i; 

G=num*f.*exp(i*num*f); 

B=fftshift(fft2(G))/(num*1 28)^2; 

IG=(B.*conj(B)); 

s(i)=sum(sum(IG(61:69.61:69)));  %  assuming  1/2  power 

contour(IG) 

%plot(s) 

%B=0:G=0;IG=0; 

end; 
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contour(IG) 

mesh(IG) 


Matlab  Function  for  calculation  N2 


h=.5:alpha=4.0948; 
lambda=476e-7; 
heff=(alpha/(1  -exp(-alpha*h))) 
wo=2. 44*  lambda*.  2/.  004; 
wo=30e-6; 

f=.406*(1-.5)^.25; 
lc=2/(pi*wo^2); 
k=2*pi/(lambda): 
l=(lc*.023)*1e-4 
shift=(1 .6)/f: 

n2=(shift/(k*l))*heff  %  as  in  Eq.  1 


Appendix  C 


Matlab  Function  for  calculation  of  Z-scan  curve 

%  March  30  1994 
hold  off, clear 
step=.5; 

load  twox  %  data  from  Z-scan  measurements  used  to  find  n2  below 
load  twoy 
Power=.001 ; 

p=max(twoy);v=min(twoy); 

PV(1)=p-v: 
n2=NTwo(Power,p,v) 
alln(1  )=n2; 

%  Curve  fitting . 

xi=.001  .*(twox(1  ):step:-twox(1 )); 
si=spline  (twox,  twoy,  xi); 
c=polyfit(twox,twoy,7); 


fit=polyval(c,twox); 

%plot(twox,twoy,xi,si,twox,fit,xi,th);grid 
%plot(twox,twoy), grid, hold  on 


for  i  =  1:2, 

r=  .00005+  i/2500  %+  i/3000 

th=DPhaser(r+.002,(p-v),xi,1);  %  DPhaser  subroutine  listed  below 

th2=DPhaser(r,(p-v),xi,0); 

plot((.001  .*twox),twoy,'g',xi,th2,'r',xi,th,'b'),hold  on 
grid 

end, pause 
hold  off 

load  threex  %  repeated  for  other  data  values 

load  threey,Power=.0025: 
p=max(threey);v=min(threey); 

PV(2)=p-v: 

n2=NTwo(Power,p,v) 

alln(2)=n2: 

si=spline(threex,threey,xi):r=. 000025 
th=DPhaser(r+.002,(p-v),xi,1):th2=:DPhaser(r,(p-v),xi,0); 
plot((.001  .*threex),threey,'g‘,xi,th2,'r',xi,th,'b’) 
grid, pause 
hold  off 

load  fourx  %  repeated  for  other  data  values 

load  foury,Power=.005; 

p=max(foury);v=min(foury);%'Fours',p,v 

PV(3)=p-v: 

n2=NTwo(Power,p,v) 

alln(3)=n2; 

si=spline(fourx,foury,xi): 

th=DPhaser(r+.002,(p-v),xi,1);th2=DPhaser(r,(p-v),xi,0); 
plot((.001  .*fourx),foury,'g',xi,th2,'r',xi,th,'b') 
grid, pause 
hold  off 

load  fivex  %  repeated  for  other  data  values 

load  fivey,Power=.01 ; 
p=max(fivey);v=min(fivey);%'Five',p,v 
PV(4)=p-v; 
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n2=NTwo(Power,p,v) 

alln(4)=n2; 

si=spline(.001  .*fivex,fivey,xi);r=.00001 5; 

th=DPhaser(r+.002,(p-v),xi,1);th2=DPhaser(r,(p-v),xi,0); 

th3=DPhaser(r,-(p-v),xi,0); 

plot((.001  .*fivex),fivey,'g',xi,th2,'r',xi,th,'r',xi,th3,'r') 
averagen2=sum(alln)/4  %  find  average  value  of  n2  from  all 

%  inputted  data 

plot(alln) 

AIW=74.8765; 

APV=sum(PV)/4; 


h=.4;alpha=6; 
lambda=51 4e-9; 
heff=(alpha/(1-exp(-alpha*h))); 
wo=2.44*lambda*.2/.004; 


f=.406*(1-.5)^.25; 


k=2*pi/(lambda*1  e2); 

shift=(APV)/f; 

n2=(shift/(k*AIW))*heff: 


Matlab  Subroutine  DPhaser 
function[T]=DPhase(radius,pv,z,a) 

c|=.44;wo=35e-6;  %  Input  beam  radius  1 .22*51 4. 5e-9*. 2/0. 0036  or 

.0015875 

k=2*pi/(514e-9);S=.005; 
if  a==0, 

a=6.4;%*(1-S)^(.35) 

else 

a=1; 

end 

zo=(k*wo^2)/2;DD=d/zo; 

Delta=pv/(0.406*(1-. 5)^(0. 25));  %  Change  sign  here  for  (de)focus 
x=(z)./zo: 

w=(wo^2)*((DD^2)*(1-(2*Delta*(DD-x).*x)./(a*DD*(1+x.^2).^2)).^2 
+(1  -(2*Delta*(DD-x))./(a*(1  +x.^2).^2)).^2);  %  as  in  Eq.  5 
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norm=(woA2)*((DDA2)  +(1  -(2*Delta*DD)./a)  .^2); 

%w=w./norm;  %  w  is  already  squared 
%  T=(1 -exp(-2*(radius'^2)./(w)))/(.5):  %  as  in  Eq.4 

TT=(1  -exp(-2*(radius^2)./(norm)))/(.5); 

T=T./mean(T); 

%hold  off, plot(x,T), pause, subplot(21 1  ),plot(x,T),grid 
%subplot(21 2), p!ot(x,T), grid,  pause, subplot(1 11) 
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applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


